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Abstract 

We have analyzed energy dissipation in a digital device ("Single-Electron 
Parametron") in which discrete degrees of freedom are used for presenting 
digital information. If the switching speed is not too high, the device may 
operate reversibly (adiabatically), and the energy dissipation £ per bit may 
be much less than the thermal energy ksT. The energy-time product £t 
is, however, much larger than Planck's constant h, at least in the standard 
"orthodox" model of single-electron tunneling, which was used in our calcu- 
lations. 
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Computation using any real physical system leads to the dissipation of energy, because 
of unavoidable coupling between the degrees of freedom which present information, and the 
environment. In most practical electronic devices (such as semiconductor transistor circuits), 
energy is dissipated at some rate even in static state, i.e. in the absence of information 
processing. Some prospective digital devices, such as various Single-Flux-Quantum0i and 
Single-Electronii logics,!! however, do not involve static power consumption, because they 
present conservative systems where digital information is coded by the choice of a local 
minimum of potential energy. In this case the energy dissipation is proportional to the 
number of logic operations. If such a conservative system is switched irreversibly (e.g., as in 
the RSFQ logic!) , the energy loss £ per one logic operation is of the order of energy barrier 
W separating the states. The barrier should be sufficiently high to make the probability p ~ 
exp(— W^/Zc^T) of thermally-induced errors low enough, so that at the physically irreversible 
computation £^min ~ ksT ^ ksT. 

Some conservative systems, e.g. the Parametric Quantron,0'i are flexible enough to allow 
independent adjustment of the potential tilt and barrier height. As has been shown earlier ,31 
this flexibility allows physically reversible (adiabatic) switching of the system with £ <^ ksT, 
if the information content of the system remains intact. If, however, information is lost 
during the operation {informationally irreversible computation), the minimum energy loss 
is at least f^min = ksT ln2 per each lost bit. Thus only at completely (physically and 
informationally) reversible computation, £ may be made much less ksT. In this case, £ 
scales as r^^, where r is the switching time, so that the product £t is fixed.lll 

£t has the dimensionality of Planck's constant, so the natural question is whether quan- 
tum mechanics imposes any fundamental lower bound on this product. A quantitative 
analysis of a particular reversible system (the Parametric Quantron) has shownii that £t 
may be much less then h. The analysis has, however, used the assumption that the potential 
energy is a function of a continuous degree of freedom (in that particular case the Josephson 
phase 0). To our knowledge, a similar analysis has never before been carried out for any 
system presenting information with discrete states. 



Recently we suggestedll2l a system ( "Single- Electron Parametron" ) based on correlated 
single-electron tunneling (see, e.g. Ref. |^), which may be used for reversible computation 
using a discrete degree of freedom - the electric charge Q. The goal of the present work was 
to find the minimum energy dissipation S for this system, and relate it to the switching time 
T and to the error probability p. We have shown, in particular, that within the "orthodox" 
theory of single-electron tunneling,0 £t is always considerably larger than h. 

Figure la shows the possible structure of a unit cell of the system. It consists of at least 
three small conducting islands (with capacitances C ^ e'^/ksT), with the middle island 
slightly shifted in y-direction. Tunnel barriers with small conductances {G <^ allow 
direct transfer of electrons only between the neighboring islands. The system is biased by a 
periodic "clock" electric field E(t) perpendicular to axis x. Let us consider the conceptually 
simplest case when the cell is charged as a whole by a single extra electron. (For practice, 
the operation using electron-hole pairs may be beneficial, but for our present discussion 
both versions are identical.) 

When the vertical component Ey{t) of the field is lower than a certain threshold value 
Et, the extra electron is kept inside the middle island. The energy diagram for this "OFF" 
state is shown at the top of Fig. lb. As Ey is increased, tunneling of the electron from the 
middle island into one of the edge islands becomes energy advantageous at Ey > Et. If the 
system is completely x— symmetric, this results in spontaneous symmetry breaking, so that 
the direction of the resulting electric dipole moment of the system is random: = Qdgff, 
dgff ^ d, Q = ±e. However, if the symmetry is broken by a weak additional external field 
E^ (say, induced by the dipole moment of a similar neighboring cell), the direction of the 
electron tunneling and hence the sign of Px is predetermined by this field. The middle frame 
in Fig. lb shows this "OFF— >0N switching" stage. Finally, when Ey is well above Ef, the 
electron is trapped inside one of the edge islands, even if the "signal" field E^ now favors 
its transfer in the opposite edge direction ("ON" state, the bottom frame in Fig. lb). In 
this state, electron transfer may only be achieved via a higher-order "co-tunneling" process; 
the probability of this process may be made negligibly small by either decreasing the tunnel 
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conductance or by inserting a few additional islands into the cell.0 If this parasitic process 
is negligible, the cell has a fixed dipole moment and may serve as a robust source of signal 
field Ex for similar neighboring cells. After this source has been used, the system is reset 
into the "OFF" state during the corresponding part of the clock cycle, when E{t) drops 
below Ef again. 

It is evident that the operation of the cell is quite similar to that of the Parametric 
Quantron,0 except now the information is presented by a discrete variable, Q = ±e. Similarly 
to the Parametric Quantron, the Single-Electron Parametron may be used for reversible 
transfer and processing of information.llilli For example. Figure 2 shows a possible structure 
of a shift register. In each neighboring cell, the extra charge sign alternates, while the 
direction of the middle island shift within plane yz is changed by O = vr — 27r/M, M > 2 
(in Fig. 2, M = 3). The clock field E()f:) has a fixed magnitude E > Et, but rotates within 
plane yz, providing periodic switching ON and OFF of the cells, with the phase shift 2tt/M 
between the neighboring cells. At an appropriate choice of E and the distance between the 
cells,0 the orientation of the dipole moment of the cells in ON state determines the direction 
of the field E^ and hence the direction of electron tunneling in the neighboring cell which 
is being switched OFF— i>ON. As a result, the information is being re-written from cell to 
cell, and thus transferred over M cells each clock period. Reversible logic operations may 
be implemented in a similar way, e.g. by using majority gates with additional output cells.§ 

Within the "orthodox" theory0 all properties of the system may be found from solution 
of the system of master equations for the probabilities Pi{t) to find the extra electron in the 
middle {i = m), left {i = I), and right {i = r) islands: 

j.Pi = {Pj^ji-Pi^ij)^ YPi{i) = 1> Pm(0) = 1, (1) 

j=m,l,r i 

where in our case the tunneling rate matrix F has only four nonvanishing components: 

, ^ ±GW , ^ ±G{W - A) 

e^l-explTW/knT]}' e^l - exp[TiW - A)/kBT]}' 

Here T^j = Tij, T^j = Tji, while W{t) ^ Ey{t)d' + const is the energy difference between 
the charge configurations with the extra electron in the middle and right islands. We will 



accept that near the decision-making point of the ON— >-OFF switching (t ^ in Fig. Ic) the 
difference is a hnear function of time: W = at. A ^ 2dEx is the energy difference between 
the left and right ON states. In order to operate with a low error probability p <^ 1, at 
the decision-making moment this difference should be large enough. Without the loss of 
generality, we may assume A > 0; then p can be found from the solution of the master 
equation as pi{oo), while the average energy S(t) dissipated by moment t can be calculated 
as 

S{t) = f_jW{t)^ + mt) - A]^}dt- (3) 

we will be mostly interested in the net dissipation S = S{oo). 

The solution of equations yields the following results. With an accuracy suf- 

ficient for our final result, the total error probability p may be calculated as a maximum 
of probabilities of the thermal and dynamic errors. The thermal error may occur due to 
thermally-activated tunneling to the wrong state (in our case, /), and its probability ptherm 
may always be expressed as exp(— A/A;bT). (Because of that, we will restrict our discussion 
to the limit A ^ ksT.) The dynamic error occurs when the switching speed a is too high, 
and the system remains in the initial (symmetric) state up to the moment when tunneling 
to the upper energy level becomes possible. If 6 = ae^ /G/^ksT ^ 1, the dynamic error 
dominates and its probability is given by expression 

p.,„ = /fo-xp(-i). K^l^- j^exp(i^)|l - Erf(^)] = 1 + 1 ^ILli^ (-,)". 

(4) 

where 7 = ae^/GA^. In order to keep pdyn <^ 1, 7 should be much smaller than 1, so that 
one can use Eq. (^ with K = 1. 

Energy dissipation depends on another dimensionless parameter, /? = ae^ /GlksT)'^ = 
(A//cbT)^7 = {A/kBT)6. Like 7 and 6, parameter j3 is also proportional to the switching 
speed a, but is much larger than both of them (because A ^ ksT) and may be comparable 
to unity. In the low-speed limit (3 <^ 1, the switching process is adiabatic. It consists of 
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numerous tunneling events (back and forth between m and r) taking place within the energy 
interval ~ ksT around the point W{t) — 0. In this case £ — K^ksT, where 



so that for this (reversible) process £ <g ksT. Notice that £ — uPhsT — nae^/GkBT 
decreases when temperature increases. 

Our model allows not only to calculate the net dissipation £, but also follow the time 
dynamics of energy transfer between the system and the environment ( "heat bath" ) during 
the switching process. During the first half of the process (when W{t) < 0) the energy 
El = — ^^(0) = T'ln2 ^ £^ is borrowed from the heat bath (which, hence, is cooled), while 
virtually the same amount S2 = S — S{0) is returned back to the heat bath during the 
second half of the process {W{t) > 0). This exchange is directly related as S{t) = TAS{t) 
to the temporal increase and consequent decrease of the entropy corresponding to the degree 
of freedom used to code information (in this particular case, the polarization P^). At the 
moment when W = 0, the system may be in either of two states (p^ = Pr = V^); i-e. 
AS — A;Bln2 has been acquired in the comparison with the definite initial state {pm — 1, 
Pr — 0). By the end of the switching (W » ksT) the entropy is restored to the initial value 
since the state is definite again (p^ = 0, pr = 1). Finite switching speed decreases Si and 
increases £2 (see the dotted lines in Fig. 3). 

In the hmit /9 ^ 1, the speed of energy change is so high that switching may take place 
only ai W > 0, but within a much larger interval of energies: AW ~ The average 

energy dissipation for this (irreversible) process is of the same order, i.e. much larger than 
ksT and independent of temperature: £ — {'K/3/2y/^kBT — (7re^Q;/2G)^/^. The results of 
numerical calculation of £ for intermediate values of /? are presented by the solid line in Fig. 



Taking into account that at p <S 1 the parameter t — A/a may be considered as the 
duration of the switching process (Fig. Ic), all our asymptotic results may be summarized 




(5) 



3. 



as follows: 
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0.67 In i , 

X < 1.97/3-1/2 




for < 1, 



for 5 < K /?, 



(6) 



2.78 (In 



2pln(l/p) 



1 



/2 , for 1< 5, /3 , 



where RQ='Kh/2e^ ~ 6.45 kf2 is the quantum unit of resistance. Since the orthodox theory 
is vahd only at GRq <^ 1, within this theory St ^ h for any switching speed. 

To summarize, we have shown that reversible computation with the energy dissipation 
S per bit well below ksT may be implemented in a physical system with discrete states. 
The quantum bound for the product St, obtained within our concrete model is, however, 
much higher than that obtained earlier for a system with continuous degrees of freedom.l^ 
Apparently the /i-limit for St may be overcome in the case of islands with discrete spectra 
of electron energies,0'0 though this may require an exponentially high energy barrier W 
during the ON state of the cell. A quantitative analysis of this opportunity is in progress. 
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FIGURES 

FIG. 1. The Single-Electron Parametron: (a) 3- island version of the system, (b) its energy 
diagram for three values of the clock field Ey, and (c) energy of the extra electron in various 
islands as a function of time, close to the decision-making moment t ^ 0. 

FIG. 2. Top (left) and side (right) views of a shift register using an array Single-Electron 
Parametron cells. Clock field E{t) rotates in yz plane. Digital bits are coded by positions of the 
extra charges in ON state of the cells, and are propagated from the top to the bottom, over M = 3 
cells during one clock period. 

FIG. 3. Components of the energy exchange between the Parametron and the heat bath as 
functions of the process speed a = dW/dt. Dotted lines: average energy flow 6i from the heat bath 
to the Parametron during the first half of the process {W < 0) and the average flow £2 from the 
device back into the heat bath during its second half {W > 0), respectively. Solid line: net energy 
dissipation £ = £2 — Si. Dashed lines show the low-speed (adiabatic) and high-speed (diabatic) 
asymptotes of the function £{a) - see formulas in the text. 
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